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Guide Field Dependence of 3D X-line Spreading During 
Collisionless Magnetic Reconnection 

L. S. Shepherd 1 and P. A. Cassak 1 

Abstract. Theoretical arguments and large-scale two-fluid simulations are used to study 
the spreading of reconnection X-lines localized in the direction of the current as a func- 
tion of the strength of the out-of-plane (guide) magnetic field. It is found that the mech- 
anism causing the spreading is different for weak and strong guide fields. In the weak 
guide field limit, spreading is due to the motion of the current carriers, as has been pre- 
viously established. However, spreading for strong guide fields is bi-directional and is due 
to the excitation of Alfven waves along the guide field. In general, we suggest that the 
X-line spreads bi-directionally with a speed governed by the faster of the two mecha- 
nisms for each direction. A prediction on the strength of the guide field at which the spread- 
ing mechanism changes is formulated and verified with three-dimensional simulations. 
Solar, magnetospheric, and laboratory applications are discussed. 



1. Introduction 

Magnetic reconnection is the basic plasma process in 
which magnetic energy is converted to kinetic and thermal 
energies (Dungey [1953]; Vasyliunas [1975]). It plays an im- 
portant role in the dynamics of explosive coronal events, 
geomagnetic substorms, solar wind coupling to the magne- 
tosphere, and magnetically confined fusion devices. Early 
models [Sweet [1958]; Parker [1957]; Petschek [1964]) and 
the predominance of numerical work on magnetic reconnec- 
tion (e.g., Birn et al. [2001]) have treated reconnection as 
two-dimensional. However, naturally occurring magnetic re- 
connection often begins in a localized region and spreads in 
the direction perpendicular to the plane of reconnection. For 
example, satellite observations of substorms in the magne- 
totail identified a dawn-dusk asymmetry caused by localized 
reconnection spreading in the westward direction (McPher- 
ron et al. [1973]; Nagai [1982]). A similar asymmetry was 
observed in the formation of arcades in the solar corona 
(Isobe et al. [2002]). Capturing effects such as these requires 
a fully three-dimensional treatment. 

A number of numerical studies have addressed X-line 
spreading in the direction of the current during quasi- 
two-dimensional reconnection. Using a magnetic per- 
turbation localized in the out-of-plane direction in Hall- 
magnetohydrodynamics (Hall-MHD), it was found that the 
localized reconnection signal propagates as a wave struc- 
ture carried by the electron current (Huba and Rudakov 
[2002, 2003]). By seeding reconnection with large random 
magnetic perturbations in Hall-MHD simulations, it was 
observed that reconnection develops into spatially isolated 
structures that lengthen in the direction of the electron cur- 
rent and that these small structures merge into larger scale 
structures (Shay et al. [2003]). It was also suggested in this 
study that spreading occurs in the direction of whichever 
species carries the current, which need not be exclusively 
electrons. Spreading by the ions when they carry the cur- 
rent was observed in hybrid simulations with localized resis- 
tivity (Karimabadi et al. [2004]). The result of these works 
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is that the reconnection X-line spreads in the out-of-plane 
direction by the current carriers in the direction of the cur- 
rent carriers (Lapenta et al. [2006]). Nakamura et al. [2012] 
presented the first systematic study to vary the fraction of 
current carried by each of the species; the results confirmed 
that X-line spreading occurs due to the current carriers. 
The results are not dependent on the Harris sheet geom- 
etry; Lukin and Linton [2011] observed X-line spreading in 
simulations of island coalescence. Note, each of these stud- 
ies primarily favored magnetotail applications, so they ei- 
ther treated anti-parallel reconnection or reconnection with 
a weak out-of-plane (guide) magnetic field compared to the 
background field. X-line spreading in a system without a 
guide field was recently observed in laboratory experiments 
at the Magnetic Reconnection experiment (MRX), and a 
physical mechanism for spreading by current carriers was 
proposed (Dorfman [2012]). 

Interestingly, experimental and satellite observations of 
systems with a strong guide field reveal strikingly different 
behavior of X-line spreading. For example, experiments per- 
formed at the Versatile Toroidal Facility (VTF) (Katz et al. 
[2010]; Egedal et al. [2011]) exhibit reconnection beginning in 
a localized region and spreading bi-directionally in the out- 
of-plane (toroidal) direction at a speed consistent with the 
Alfven speed based on the guide field. Another example is 
bi-directional spreading (or elongating) of ribbons observed 
during two-ribbon solar flares (Qiu [2009]), including the 
Bastille Day flare (Qiu et al. [2010]). This presumably is 
related to spreading of the looptop reconnection site where 
a sizable guide field is likely to be present. This spread- 
ing was also inferred to take place at the local Alfven speed. 
Prominence eruptions in the corona have also been observed 
to spread bi-directionally; this behavior was attributed to 
magnetic reconnection propagating along the magnetic po- 
larity inversion line (PIL) (Tripathi et al. [2006]). In magne- 
tospheric contexts, observations of extended X-lines several 
Earth radii long at the magnetopause (Phan et al. [2000]; 
Fuselier et al. [2002]) and hundreds of Earth radii in the 
solar wind (Phan et al. [2006]) suggest that X-line spread- 
ing occurs in these areas as well, although direct evidence of 
spreading is prohibitively difficult with single- or even multi- 
point satellite observations. X-line spreading was also seen 
in three-dimensional two-fluid simulations with a guide field 
(Schreier et al. [2011]). 

The existing observational data provide a clear indication 
that the mechanism controlling X-line spreading strongly de- 
pends on the strength of the guide field. In the weak guide 
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field limit, the signal is transmitted by the current carri- 
ers; in the strong guide field limit, the reconnection signal 
is transmitted by the magnetic field as an Alfven wave. We 
hypothesize that, in general, the X-line spreads in both di- 
rections at the speed of whichever mechanism is faster for 
that direction. In this paper, we present an estimate of the 
critical guide field where the spreading mechanism changes 
and confirm the theory with three-dimensional two-fluid nu- 
merical simulations. 

The layout of this paper is as follows. A prediction of 
the critical guide field at which the spreading mechanism 
changes from current carriers to Alfven waves is developed 
in Sec. 2. The simulation setup and results are discussed in 
Sees. 1 and 4, respectively. A discussion of the results and 
potential applications is in Sec. 5. We emphasize that we are 
considering current sheets that are already thin, with large 
amounts of free magnetic energy present. The important 
topics of how the sheets become thin and how the magnetic 
energy is stored is outside the scope of this paper. 

2. Theory 

Here, we develop a prediction of the speed at which the 
X-line spreads in each out-of-plane direction as a function 
of guide field and derive the critical guide field at which 
the mechanism causing the spreading changes from current 
carriers to Alfven waves. To do so, we make the following 
simplifying assumptions. We treat a quasi- two-dimensional 
system, meaning that the equilibrium parameters do not de- 
pend strongly on the direction normal to the reconnection 
plane for all time. We assume the current layer is flat, so 
that the current sheet is either not curved or that the cur- 
vature does not strongly contribute to the dynamics. We 
assume the plasma parameters are symmetric on either side 
of the current layer; asymmetries (Cassak and Shay [2007]) 
are not considered here. Finally, we assume that a single 
mode dominates the dynamics; in previous simulations, it 
was shown that when multiple modes of reconnection occur, 
they can impede the spreading of X-lines (Schreier et al. 
[2011]). This assumption is valid at early times and in sys- 
tems in which only a single mode is present. 

First, we estimate the spreading speed in each direction 
for each spreading mechanism. We begin with the speed due 
to the current carriers. From Ampere's law, the current is 
J = cV x B/47T, where B is the magnetic field. For sim- 
plicity, we first assume the electrons carry the out-of-plane 
current, so that the electron velocity is v e = — J/ne, where 
n is the electron density and e is the proton charge. Using a 
scaling argument, the electron speed v eg in the out-of-plane 
direction is 
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where B rec is the strength of the reconnecting magnetic field 
upstream of the electron layer, 5 is the thickness of the cur- 
rent layer, and g refers to the direction of the guide field. As 
has been previously established (Huba and Rudakov [2002]; 
Shay et al. [2003]; Karimabadi et al [2004]; Lapenta et al. 
[2006]; Lukin and Linton [2011]; Nakamura et al. [2012]), 
this is the X-line spreading speed in the absence of a guide 
field. In the strong guide field limit, the observations sug- 
gest the spreading speed is the Alfven speed CA g based on 
the guide field, given by 



CAg 



B a 



7T 771*71 



(2) 



where B g is the strength of the guide field and vra is the 
proton mass. 

Our hypothesis is that the X-line spreading speed in the 
direction of the electron out-of-plane flow, which we call vx e . 



is the larger of v eg and CA g ' 

v Xe - max{7J e p,CAs}. 



(3) 



From this, one can find the critical guide field B cr it, e at 
which the spreading mechanism changes, where the e sub- 
script denotes the critical field for motion in the direction 
of the out-of-plane electron flow. Setting Eq. (1) equal to 
Eq. (2) and solving for B g gives 
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where di = c/uo V i is the ion inertial length and uo P i — 
(Aime 2 /mi) 1 / 2 is the ion plasma frequency. Since S is typ- 
ically less than di as the current is set by electron scales, 
we expect B cr it,e > B rec , although B cr it,e is on the order 
of and slightly larger than the reconnecting magnetic field 
strength upstream of the ion dissipation region. 

We perform a similar analysis for the spreading speed 
in the direction of the ions vxi- Since electrons carry the 
current, the ion speed v% 9 in the out-of-plane direction is 
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Figure 1. Schematic diagram showing the mechanisms 
that cause X-line spreading. The thin blue arrows are 
the reconnecting magnetic field components, the yellow 
arrow is the total current. The red arrows denote the 
speed of the current carriers; the thick blue arrows de- 
note the speed of Alfven waves along the guide field. The 
top, middle, and bottom diagrams show the spreading 
mechanisms for strong, weak, and arbitrary guide field 
strengths. In each case, X-line spreading occurs at the 
faster speed in each direction. 
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Therefore, the X-line spreading speed in the direction of the 
ion out-of-plane flow vx% — max{vi g , CA g } is given by the 
Alfven speed based on the guide field, 
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(6) 



Since vi g — 0, the critical guide field B cr it,i for spreading in 
the direction of the ion out-of-plane flow is 



B C rit,i — 0. 



(7) 



These results can be generalized to systems with both 
electrons and ions carrying some of the current. Following 
Nakamura et al. [2012], we define the fraction of the total 
current J z carried by the ions as a, which is assumed known 
or measurable. Letting J iz = aJ z , one has J ez = (1 — a) J z 
so that J z = Ji Z + J ez - By performing a similar analysis 
as before, one finds the out-of-plane electron and ion flow 
speeds due to current carrying are 



~ (1 — a) 
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which generalizes Eqs. (1) and (5). The X-line spreading 
speeds in the direction of the electron and ion out-of-plane 
flow are 



V X e = max{v eg , CAg}, 

v X i - max{v ig ,CAg}, 



(9) 



respectively, which generalizes Eqs. (3) and (6). Finally, 
the critical guide fields at which the mechanism for X-line 
spreading changes from the current carriers to Alfven waves 
in the direction of electron and ion flows are given by 



B C rit,e ~ (1 — a)B rec -r, 




B C rit,i ^ OtB rec e , 




(10) 



respectively, which generalizes Eqs. (4) and (7). 

The predictions derived here are summarized pictorially 
in Fig. 1, where the current is depicted by the yellow arrows 
and the reconnecting magnetic fields are the thin blue lines. 
The thick arrows denote the speeds of the current carriers 
(in red) and the Alfven speed (in blue) in each out-of-plane 
direction. The top, middle, and bottom plots show the re- 
sults for strong, weak, and arbitrary guide field strengths, 
respectively. In each case, the X-line spreading speed is the 
longer of the arrows on either side. We point out that there 
is nothing preventing the spreading mechanisms from be- 
ing different in the two directions, i.e., Alfven waves in one 
direction and current carriers in the other, if that is what 
Eq. (10) dictates for the system parameters. 

3. Simulation Setup 

To test the predictions on X-line spreading, three- 
dimensional numerical simulations are performed using the 
two-fluid code F3D (Shay et al. [2004]). The code updates 
the continuity, momentum, and induction equations with 
the generalized Ohm's law including electron inertia. Mag- 
netic fields and densities are normalized to arbitrary values 
Bo and no. Velocities are normalized to the Alfven speed 
cao — Bo I (4n rriino) 1/ ' 2 . Lengths are normalized to the ion 
inertial length dio = c/u P io — (rriiC 2 /Annoe 2 ) 1 ^ 2 . Times are 



electric fields to Eo = caoBo/c, and temperatures to To = 
rrnc 2 Ao. 

Simulations are performed in a three-dimensional domain 
of size L x x L y x L z = 51.2 x 25.6 x 256.0 dm, where x is the 
direction of the oppositely directed field, y corresponds to 
the inflow direction if the simulations were two-dimensional, 
and z is the direction of the initial current. The plasma is 
assumed to be isothermal and there is no resistivity (77 = 0). 
Boundaries in all three directions are periodic, but the sys- 
tem is long enough in the z direction that the periodic 
boundaries do not affect the dynamics on the time scales 
of import to the present study. 

For simplicity, the simulations have the electrons carrying 
all of the initial current (i.e., a — 0). The electron inertia 
is m e — mi/25. In previous simulations with this electron 
mass (and confirmed in the simulations here), it has been 
observed that the current layer thickness 5 thins down to 
the electron inertial scale d e — 0.2 di and the reconnecting 
magnetic field at the electron layer is B rec — 0.4 Bo (Jemella 
et al. [2003]). Substituting this into Eq. (4), we predict a 
critical guide field of 



: 2Bo 
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normalized to the ion cyclotron time Q c 
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Therefore, we can test the theory by running a series of 
simulations in which the initial guide fields are B g = 
0,0.5,1,1.5,2,2.5 and 3. Note, the scaling 5 ~ d e and 
B rec ~ 0.4 Bo may or may not be representative of naturally 
occurring reconnection; care should be taken to investigate 
this for particular applications. 

The initial configuration is a double tearing mode with 
two Harris sheets, B x o(y) = taxLh[(zj + L y /4) / wo] — tanh[(y — 
L y /4)/wo] — l, with uniform initial temperature T — 1 and a 
non- uniform plasma density to balance total pressure. Here, 
wo = 0.4 dio is the initial current layer thickness. We choose 
this thickness to be comparable to the smallest value of the 
ion Larmor radius p s — c s /Q C i — VT/B g ~ 0.33, where c s 
is the sound speed and the latter expression is written in 
normalized units. This scale is the Hall scale in the pres- 
ence of a strong guide field (Zakharov et al. [1993]; Rogers 
et al. [2001]). It is worth noting that the Hall scale increases 
smoothly from p s to d% as the guide field is decreased to 
zero, which follows from a linear analysis of Hall-MHD waves 
(Rogers et al. [2001]). Consequently, the smaller guide field 
simulations start with a current sheet that is thin relative 
to the Hall scale, and should onset rapidly. As the guide 
field is increased, the time to onset should increase and it 
is expected that a hyper-resistive phase of reconnection will 
occur before onset. This behavior will not adversely impact 
our study, as we will separate out the times for which Hall 
reconnection is dominant. 

We employ a grid scale of Ax x Ay x Az = 0.05 x 
0.05 x 1.0 dio. Using a stretched grid in the out-of-plane 
direction has been done before (Shay et al. [2003]), and is 
acceptable since the in-plane kinetic-scale dynamics is on 
smaller scales than the out-of-plane dynamics. To ensure 
the stretched grid scale in the out-of-plane direction does 
not play a role in the numerics, some simulations are con- 
firmed by comparison with simulations with Az = 0.5 dio. 
All equations employ a fourth-order diffusion with coeffi- 
cient Dax — Da v — 2.5 x 10 ~ 5 in the x and y directions. 
In the out-of-plane direction the fourth-order diffusion co- 
efficient D^g depends on the speeds in the out-of-plane di- 
rection. For B g < 2.0 the fourth-order diffusion coefficient 
is D^g — 0.064 and for B g = 2.5 and 3.0 the fourth-order 
diffusion coefficient is D± g — 0.081 and 0.097, respectively. 
The values of D± g were tested by varying the value by a 
factor of two to ensure that D^ g does not play a significant 
role in the dynamics. 

The inclusion of a guide field in these simulations changes 
the nature of reconnection relative to previous work on X- 
line spreading. In three-dimensional periodic domains, it is 



X - 4 



SHEPHERD AND CASSAK: X-LINE SPREADING 




-100 -50 




z 



50 100 



Figure 2. The envelope f z (z) used to localize the mag- 
netic perturbation in the out-of-plane direction z. 



well established that the linear tearing instability is excited 
where k • B = 0, where B = B x0 St + B z oz is the equilib- 
rium magnetic field and k = k x St-\- k z z is the wave vector of 
the mode. The periodic domain enforces that k x = 2nm/L x 
and k z = 2nn/ 'L z , where m and n are integer mode numbers 
which specify the number of X-lines in the x and z direction, 
respectively. In the absence of a guide field, this condition 
is only satisfied where B x o = 0. With a guide field, it is 
satisfied wherever q(y) = L x B z (y) / L z B x (y) = m/n is a ra- 
tional number, where q(y) is the safety factor well known 
in fusion applications. The y locations where k • Bo = is 
satisfied are called rational surfaces, and for the equilibrium 
profile here, the modes are displaced from where B x0 = by 
a distance y s — wotanh -1 (nL x B g /mL z Bo). Thus, modes 
in our simulations are excited on multiple rational surfaces. 

Reconnection is seeded using a coherent magnetic pertur- 
bation localized in the out-of-plane direction of the form 



4. Results 

To ensure numerical feasibility of the simulations in three 
dimensions, we benchmark the simulations in two dimen- 
sions. The simulations are evolved from t — until non- 
linear reconnection develops. In two dimensions, symmetry 
dictates that the n — mode is the only excited mode. The 
reconnection rate E, measured as the time rate of change of 
the difference in magnetic flux between the X-line and O-line 
during a quasi-steady period, is approximately 0.08-0.1 for 
all simulations. Also, as expected, the time until Hall recon- 
nection begins increases as the guide field increases since wo 
is held fixed, and there is a brief hyper- resistive reconnection 
phase before onset for stronger guide fields. 

In the three-dimensional simulations, the evolution at 
z — is very similar to what is observed in two dimen- 
sions: the time scale of the development of reconnection is 
comparable, and a hyper-resistive phase precedes Hall recon- 
nection. The reconnection rates can be compared, as well. 
The reconnection rate in three-dimensions is measured by 
taking a cut of the z component of v x B in the ^-direction 
across the X-line; far from the current sheet, it asymptotes 
to the reconnection electric field in a steady-state. The re- 
connection rates are in the 0.08-0.1 range, comparable to the 
two-dimensional results. One noticeable difference, as dis- 
cussed earlier, is that oblique modes dominate over n = 
modes in three-dimensions. 

Each of the three-dimensional simulations display some 
form of X-line spreading. This can be seen qualitatively in 
Fig. 3 for the simulation with B g = 3. The out-of-plane cur- 
rent J z at time t = 30 (top row) and at time t = 35 (bottom 
row) is displayed at three different out-of-plane positions: 
z = —30, 0, and 30 from left to right. At the earlier time 
t = 30, a transition to fast (Hall) reconnection at multiple 
sites at z = has occurred, consistent with the development 
of multiple oblique modes. At z = ±30, the reconnection 
is still hyper-resistive. At the later time t = 35, the cur- 
rent sheet at all three positions in z has developed multiple 
Hall reconnection X-lines. Thus, the Hall reconnection sig- 



B yl = B 1 sin(k x x + k z z)f z (z), (12) nal propagates bi-directionally from z = for B g = 3. It 



where B\ — 0.1. Here, f z (z) is an envelope that localizes 
the perturbation in the out-of-plane direction and is given 
by f z (z) {tanh[(z + w 0z )/6] - tanh[(z - w 0z )/6]}/2. We 
use wo z = 1; a plot of f z (z) is in Fig. 2. Random magnetic 
perturbations that range from m,n — to 20 with small 
amplitude 0.02 Bo are included with the initial conditions 
to break symmetry so that secondary islands are ejected. 

In early simulations, we initially excite only the (m, n) = 
(1,0) mode in Eq. (12). Even though this mode is the 
strongest perturbation, oblique modes with n/0 grow from 
the noise and dominate the reconnection. This is consis- 
tent with recent particle-in-cell (PIC) simulations (Daughton 
et al. [2011]) and linear theory (Baalrud et al. [2012]). 
Oblique modes in reconnection have been observed many 
times in fusion applications (see e.g., Grasso et al. [2007]). 
In light of these results, we include oblique modes in Eq. (12) 
and compare the results with the original simulations. The 
values of m and n are chosen so that the displacement y s 
is less than wo. In this study, m — 1 for all simulations 
and n ranges from to 3. Initially exciting oblique modes 
has no noticeable effect on the results on the development 
of reconnection. Thus, the results of this study are expected 
to be independent of the modes used to seed reconnection. 
Note, although the modes are oblique, they are still quasi 
two-dimensional until they start interacting strongly. It was 
shown (Schreier et al. [2011]) that interacting oblique modes 
can prevent X-lines from spreading, so we focus on times 
early enough in the evolution that oblique mode interactions 
have not yet occurred. 



is worth noting that the multiple oblique mode reconnec- 
tion seen here is consistent with previous simulations, and 
the reason multiple X-lines appear despite the m — 1 mode 
being the dominant mode is that there are multiple modes 
simultaneously excited on different rational surfaces. 

To quantify the speed at which the X-line spreads, we 
must develop a systematic way to determine the extent of 
the reconnection region. As Hall reconnection develops, the 
out-of-plane current J z at the X-line becomes noticeably 
higher than regions where reconnection is hyper-resistive. 
For each slice in z, we measure the maximum out-of-plane 
current, which we call Jmax(z). These maximum values of 
the current correspond to the location of the X-line for each 
position in z. The extent of the X-line can then be readily 
seen in a stack plot of Jmax(z) as a function of t. Stack plots 
for all six initial guide fields in this study are displayed in 
Fig. 4. As mentioned earlier, the plots only cover early times 
when the three-dimensional X-line structure is well defined 
because the interaction of oblique modes make defining the 
X-line structure prohibitive. 

The bright white regions in Fig. 4 correspond to the 
strongest currents and, thus, the Hall reconnection X-lines. 
The dimmer areas outside of the white dashed lines (the 
red to black colors) indicate the region undergoing hyper- 
resistive reconnection. As expected, the B g — simulation 
onsets almost immediately without a hyper-resistive phase 
since wo < di, while the onset time increases as the guide 
field increases, leading to a longer hyper-resistive phase. 
Both phases of reconnection spread in the z direction as 
time evolves; we focus on the Hall reconnection X-lines in 
the present study. 
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Figure 4. Stack plots of J ma x(z) as a function of t and z. The vertical dashed green lines indicate the 
range of time over which the spreading speed is measured, U and tf. The dashed white lines indicate the 
extent in z of the X-line; their slope gives the speed of the spreading. Note, the images in the bottom 
row are on a different scale in z than those in the top row. 



From Fig. 4, the qualitative differences in the nature of 
the spreading as a function of the guide field can readily be 
seen. For the strong guide field simulations B g > 2, the 
X-line spreads symmetrically about z = (top row), which 
is consistent with our expectations for the strong guide field 
regime from Eq. (11). However, for simulations with a guide 
field weaker than the predicted condition B g <2,we observe 
different spreading behavior in the +z and — z directions 
(bottom row). For the B g = 1.5 case, there is bi-directional 
spreading, as observed in the stronger guide field runs, but 
the spreading is not symmetric about z = 0. The spreading 
in the — z direction appears marginally faster than in the 
+z direction. These differences are further amplified in the 
B g = 1 simulation. With no guide field (B g = 0), spreading 
occurs primarily in the — z direction, with negligible spread- 



ing in the +z direction. Since J z is in the -\-z direction for 
this reconnection site, the propagation is in the direction of 
the electron out-of-plane flow, consistent with previous work 
(Huba and Rudakov [2002]). 

To make this quantitative, we measure the spreading 
speed of the X-line after Hall reconnection begins by find- 
ing the length of the X-line in the out-of-plane direction; 
its time rate of change between an initial and final time is 
the spreading speed. To do so, we note that the reconnec- 
tion rate during hyper-resistive reconnection never exceeds 
0.01. We observe that when the out-of-plane Hall electric 
field Eu g = (J x B) p /nec in a cut in the y direction through 
the X-line exceeds 0.01, the reconnection has begun its tran- 
sition to Hall reconnection. We also note empirically that 
Jmax at the time of this transition is always close to 3.3, 
which is robust for all the simulations performed here. Thus, 
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Figure 5. Cuts of the stack plot for the B g = 2.5 case 
plotted in Fig. 4. The dashed (red) line and solid (blue) 
line are at U = 18 and tf = 25, respectively. The hor- 
izontal dotted line is at J thresh as defined in the text. 
The vertical dotted lines mark z = ±5, the approximate 
extent of the initial magnetic perturbation. The green 
line denotes the change in length of the X-line between 
ti and tf. 



we take Hall reconnection as occurring when Jmax exceeds 
a threshold value of Jthresh = 3.3. 

The time frame over which the spreading speed is mea- 
sured is defined as follows. The initial time U is defined as 
the earliest time that Jmax exceeds Jthresh over the entire 
range from z = d=5. This range of z is chosen because the 
initial magnetic perturbation that seeds the X-lines is local- 
ized in this region, so genuine spreading not being influenced 
by the growth of reconnection inside the initially perturbed 
region requires the signal to leave this range in z. The final 
time t/ is defined for each simulation as the latest time in 
the evolution before multiple oblique modes interact; this 
assessment is done visually by finding where the current de- 
velops complicated structure as seen in Fig. 3. The length 
of the X-line at a given time is defined as the extent in z for 
which Jmax exceeds Jthresh- The spreading speed is calcu- 
lated as the difference of the length of the X-line between tf 
and ti divided by the time difference. 

An example of this procedure is presented in Fig. 5, where 
representative data for the B g — 2.5 simulation is shown. 
The initial time is ti = 18, which is the earliest time that 
Jmax > Jthresh everywhere between z — ±5, as shown by 
the red dashed line. The final time is taken to be tf — 25; 
a plot of Jmax(z) at t f is shown as the blue line. The hori- 
zontal dotted line marks the current threshold Jthresh = 3.3 
and the vertical dotted lines mark the boundary of z = ±5. 
The change in length between the two times is the distance 
between the curves at Jthresh, marked by the green line seg- 
ments. The lengths and speeds are calculated separately for 
the ±z direction because the speeds in the two directions 
may be different depending on the strength of the guide 
field. For the B g — 2.5 simulation, the change in length in 
the +z and — z directions are 21 and 20 dio, respectively, and 
dividing by the time difference gives speeds of vxi = 3.0 cao 
for the speed in the +z direction (the direction of ion out- 
of-plane flow) and vx e = 2.9 cao for the speed in the — z 
direction (the direction of electron out-of-plane flow). 

The initial and final times U and tf for each simula- 
tion are illustrated in Fig. 4 as the vertical dotted green 



lines. The dashed white lines connect the extent of the X- 
line at the initial and final times. By inspection, one can 
see that the technique we employ to measure the extent of 
the X-line appropriately captures the evolution of the X-line 
length. Also, since the region of stronger current is rather 
straight between the beginning and final times, this implies 
the spreading speed is approximately constant in time. 

The measured X-line spreading speeds vx e and vxi are 
calculated as the time rate of change of the length of the X- 
line, which is equivalent to the slope of the white dashed lines 
in Fig. 4. The results for the spreading speed in both direc- 
tions are plotted as a function of guide field B g in Fig. 6. 
The measured value of the spreading speed is given by the 
solid blue triangles for vxi and the red stars for vxe- Note, 
vxi for B g = 0.0 and 0.5 is plotted as zero as the hollow 
blue triangles. This is because the Hall reconnection signal 
is found to not extend past z = ±5 for either simulation 
during the time considered. 

To compare these results to the theory, note that the 
electrons carry all of the out-of-plane current in the — z di- 
rection in our simulations. Therefore, in the weak guide field 
regime B g < 2, Eq. (3) predicts that the spreading speed 
in the direction of the electron current vxe is the speed of 
the electrons given in Eq. (1), which is independent of B g . 
When B g > 2, the spreading speed is determined by the 
Alfven speed given by Eq. (2), which increases linearly with 
B g . The predicted speed of X-line spreading in the direction 
of the ion current vx% is the Alfven speed due to the guide 
field, as given by Eq. (6), which increases linearly with B g 
for all guide field strengths. 

The predicted spreading speeds vx e and vxi are depicted 
in Fig. 6 by the solid red line and the dashed blue line, 
respectively. Qualitatively, the data reveal that the nature 
of X-line spreading is sensitive to the strength of the guide 
field. To interpret this more quantitatively, we first discuss 
the estimated uncertainties in our speed measurements. If 
we use a higher value of the current threshold Jthresh, the 
spreading speed changes on the order of 15-20%, which we 
take as the uncertainty. We note that for the large guide 
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Figure 6. Spreading speeds vx% and vxe as a function 
of guide field B g . The red asterisks are the measured 
values of vxe] the solid red line is the prediction from 
Eq. (3). The solid blue triangles are the measured values 
of vx%] the open triangles are for simulations for which 
no spreading was measured. The dashed blue line is the 
prediction from Eq. (6). 
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field runs B g > 2, the speeds in either direction are within 
the uncertainties of each other. However, for B g < 2, the 
speeds in either direction are separated by more than their 
uncertainty. These two results suggest that the spreading 
mechanism is the same in both directions for B g > 2 and is 
different in either direction for B g < 2, which quantitatively 
agrees with Eq. (11). 

For the absolute spreading speeds, when the estimated 
uncertainties are taken into account, the measured values 
agree pretty well with the predicted speeds. It is unexpected 
that the B g = 3 speeds are slower than B g = 2.5, but both 
are within the uncertainties of the predicted value. Also, 
it is expected that for B g = 0.5, a non-zero value could be 
obtained if there had been a longer time before the oblique 
modes started interacting. Therefore, we conclude that data 
in Fig. 6 quantitatively support the theory presented in 
Sec. 2. 

In conclusion, the mechanism of X-line spreading in the 
out-of-plane direction is qualitatively different depending on 
the strength of the guide magnetic field. For B g > B cr it,e, 
X-line spreading occurs bi-directionally along the guide field 
at the Alfven speed. For B g < B cr i t , e , X-line spreading oc- 
curs bi-directionally along the guide field, but the spreading 
speed in the direction of the current carriers is the speed 
of the current carriers and in the direction opposite of the 
primary current carriers the spreading speed is the Alfven 
speed. Measurements of X-line spreading for the hyper- 
resistive reconnection that precedes Hall reconnection agree 
with the results obtained from measuring the Hall reconnec- 
tion spreading (not shown). Therefore, the main result of 
this study applies both to Hall and hyper- resistive reconnec- 
tion in a two-fluid model. 

5. Discussion 

In summary, the mechanism of X-line spreading in the 
out-of-plane direction is qualitatively different for strong 
guide magnetic fields than it is for weak guide fields. For 
weak guide fields, the reconnection signal is propagated by 
the current carriers, as has previously been established; for 
strong guide fields, the reconnection signal is propagated by 
Alfven waves along the guide field. In general, the spreading 
speed in either out-of-plane direction is given by the max- 
imum of the speed of the current carriers in that direction 
and the Alfven speed based on the guide field, as given by 
Eq. (9). 

Because the changeover from one spreading mechanism 
to the other is abrupt, there is a critical guide field strength 
(for each direction) at which the nature of the spreading 
switches. This critical field depends only on the strength of 
the reconnecting magnetic field, the ion inert ial scale, the 
thickness of the electron dissipation region, and the frac- 
tion of the current carried by each species, as given by 
Eq. (10). When the guide field B g exceeds the critical field, 
the spreading is due to Alfven waves; when it is smaller, the 
spreading is due to the current carriers. The weak guide 
field result is consistent with previous numerical work of X- 
line spreading (Huba and Rudakov [2002]; Shay et al. [2003]; 
Karimabadi et al. [2004]; Lapenta et al [2006]; Nakamura 
et al. [2012]), but the new result generalizes the predictions 
to include a guide field. 

The present results may be relevant for interpreting ob- 
servations of reconnection in many settings. For example, in 
laboratory experiments, X-line spreading has been observed 
to be bi-directional and at the Alfven speed in the strong 
guide field limit (Katz et al. [2010]) and uni-directional in 
the small guide field limit (Dorfman [2012]). These results 
are consistent with the results of the present study. 

Another potential application is for solar flares. Two- 
ribbon flare evolution is marked by the ribbons moving 
apart from each other as time evolves, which is interpreted 



as newly reconnected field lines piling on top of previously 
reconnected field lines. In addition to this behavior, bi- 
directional spreading or elongation of the ribbon in the di- 
rection parallel to the ribbons along the polarity inversion 
line has been observed (Qiu [2009]). It was shown that the 
spreading speed was consistent with the Alfven speed (Qiu 
[2009]). Since the reconnection driving the flare is most 
likely to have a sizable guide field, the present results sug- 
gest that this type of bi-directional spreading at the Alfven 
speed would be expected. 

The Bastille Day flare exhibits this spreading, as well 
(Qiu et al. [2010]). From geometrical considerations of the 
magnetic fields of the flare loops, it was argued that the 
guide field was of comparable size as the reconnecting field, 
with B g ~ 0.4 — 1.2 times the reconnecting field (Qiu et al. 
[2010]). We can check this using the present results and the 
observed properties of the spreading. From the observations, 
the spreading speed ranged between 30—70 km/s (Qiu et al. 
[2010]). Let us assume the spreading is governed by Alfven 
waves. Assuming an average density of n = 10 13 cm -3 (Qiu 
et al. [2010]), the guide field ranges from B g ~ 15 - 100 G 
using Eq. (2). The motion of the ribbons normal to the 
ribbons was 20 km/s (Qiu et al. [2010]), which is expected 
to be correlated to the inflow speed at the reconnection site. 
Since the inflow speed is often taken to be 0.1 of the Alfven 
speed based on the reconnecting field, the reconnecting field 
strength B rec ~ 140 G. These results suggest the guide field 
is about 0.1-0.7 of the reconnecting field. Despite the large 
uncertainties, the two techniques give similar results. This 
analysis is obviously oversimplified and merely presented as 
an example of how the results can be used, but it is hoped 
that future work will allow for a meaningful assessment of 
the relative strengths of the guide and reconnecting fields. 
The reason this may be useful, as emphasized by Qiu et al. 
[2010], is that the strength of the guide field is known to 
influence the production of secondary islands (Drake et al. 
[2006b]), and it has been suggested that the presence of 
secondary islands (plasmoids) is important for particle ac- 
celeration (Drake et al. [2006a]). 

Another interesting application is for reconnection in the 
solar wind, where reconnection X-lines extending hundreds 
of Earth radii have been reported (Phan et al. [2006]). One 
can ask whether the spreading of reconnection in the out- 
of-plane direction could allow the X-line to be that long. To 
estimate the size of X-lines, assume that reconnection begins 
close to the Sun with an initially small finite length in the 
out-of-plane direction. Suppose the reconnection site con- 
vects out with the solar wind at a speed vsw- (For simplic- 
ity, this calculation ignores variations in solar wind speed, 
magnetic field strength, and plasma density as a function of 
distance from the Sun.) Then, the time it takes to get to a 
position 77 away from the Sun is t ~ rf/vsw- If the speed 
of the spreading of the X-line is vx, then the extent L of 
the X-line at 77 is L ~ vxt ~ rfVx /vsw, which gives the 
upper limit on the length of the X-line that could arise in 
the solar wind. 

One can test the implications of this from the observa- 
tions of the reconnection event in the Phan et al. [2006] 
study, where a solar wind speed is inferred to be vsw = 
340 km/s. The satellite observations occurred near the 
Earth, so 77 ~ 1 AU ~ 2.3 x l^R E - The Alfven speed 
based on a guide field of strength B g — 4 nT and density 
n = 20 cm 3 (Phan et al. [2006]) is c Ag = 19 km/s. If we take 
this as the spreading speed vx, then the maximum length 
of the X-line is L ~ rfVx/vsw ~ 1.3 x 10 3 Re- This exceeds 
the length of the X-line reported by Phan et al. [2006] , which 
was 390 Re- In this event, the strength of the guide field 
was 0.35 of the reconnecting field, so the Alfven speed is the 
slower of the two velocities and the extent of the X-line if 
spreading is due to current carriers is even longer. Thus, 
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while this calculation assumes that the reconnection pro- 
ceeds at short distance from the Sun, it gives an indication 
that it is not impossible to achieve reconnection X-lines of 
the lengths reported by Phan et al [2006]. Further work is 
necessary to make more careful comparisons of the theory 
to data. 

We conclude by collecting some assumptions of this work. 
We treat our system as quasi- two- dimensional, meaning any 
variation in the system in the direction of the current is 
negligible. The current sheet in all simulations performed 
are initially thin, meaning that free magnetic energy has 
already been stored. The plasma parameters across the cur- 
rent sheet are assumed symmetric. The simulations employ 
a two-fluid model, which does not fully capture electron scale 
physics. This may make quantitative changes to our results 
(such as the estimate of 6 and the size of B rec ), but we 
do not expect qualitative changes to the theoretical results. 
Also, the simulations are isothermal and contain no thermal 
conduction. 

Another main assumption is that only a single mode is 
dominating the dynamics. However, the role of multiple 
oblique modes can play an important part of the dynam- 
ics of the spreading process. As seen in our simulations, 
the X-line structure is identifiable at early times but as the 
complicated nature of the oblique modes develop, the X-line 
structure break up due to the interaction between the cur- 
rent sheet. The interaction of oblique modes can impede 
X-line spreading (Schreier et al. [2011]). More work is nec- 
essary on this topic. 
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